Avalanche Dynamics in Wet Granular Materials 
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A detailed characterization of avalanche dynamics of wet granular media in a rotating drum 
apparatus is presented. The results confirm the existence of the three wetness regimes observed 
previously: the granular, the correlated and the viscoplastic regime. These regimes show qualita- 
tively different dynamic behaviors which are reflected in all the investigated quantities. We discuss 
the effect of interstitial liquid on the characteristic angles of the material and on the avalanche size 
distribution. These data also reveal logarithmic aging and allow us to map out the phase diagram of 
the dynamical behavior as a function of liquid content and flow rate. Via quantitative measurements 
of the flow velocity and the granular flux during avalanches, we characterize novel avalanche types 
unique to wet media. We also explore the details of viscoplastic flow (observed at the highest liquid 
contents) in which there are lasting contacts during flow, leading to coherence across the entire 
sample. This coherence leads to a velocity independent flow depth at high rotation rates and novel 
robust pattern formation in the granular surface. 

PACS numbers: 45.70.-n, 45.70.Ht, 45.70.Mg 
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INTRODUCTION 



While most research on the physics of granular media 
has focused on dry grains, the presence of even micro- 
scopic quantities of interstitial liquid can have profound 
effects on the physical behavior |, ||, |, | @, |, |. 
Capillary forces lead to cohesion which greatly enhances 
the stability of wet samples, and several recent stud- 
ies have investigated liquid induced effects on the static 
properties [|[ |[ |[ |[ |[ 0| . In previous studies of the re- 
pose angle of wet granular media using a draining crater 
apparatus JlOt , we identified three fundamental regimes 
as a function of the liquid content. The granular regime 
at very low liquid contents is dominated by the motion of 
individual grains; in the correlated regime corresponding 
to intermediate liquid contents, a rough surface is formed 
by the flow of separated clumps; and the repose angle of 
very wet samples results from smooth, cohesive flow with 
viscoplastic properties. The addition of liquid also quali- 
tatively changes the dynamic behavior of granular media, 
as was observed qualitatively in our angle of repose mea- 
surements and also revealed in recent segregation studies 
^. A quantitative characterization of the dynamics 
of wet granular flow, however, has thus far been miss- 
ing. In this paper we study the avalanche dynamics and 
flow properties of wet granular materials to investigate 
how the gradual addition of a small amount of liquid in- 
fluences the dynamic properties of granular surface flow, 
and we find that the three fundamental regimes observed 
earlier are also reflected in the present experiment. Some 
of the results presented here have been published else- 
where [ pT| . 

Granular surface flow occurs when the inclined surface 
of a granular medium loses its stability against gravita- 



tional force. There are several experimental methods to 
investigate such flow, e.g. by gradually tilting a granular 
sample, which leads to an avalanche when the surface an- 
gle exceeds a critical value 9^ax IHl ■ Alternatively, if the 
surface angle is close to 0max , avalanches can be triggered 
by point-like perturbations The inclination of the 

surface can also be increased by the addition of grains to 
the top of the pile |14, ^ or by slowly decreasing the 
size of the supporting base . 

In our measurements we employ a rotating drum appa- 
ratus (a cylindrical chamber partly filled with a granular 
medium and rotated around a horizontal axis) | pT| , p^ , 
which tilts the sample in a highly controlled and repro- 
ducible manner and offers several benefits. The most im- 
portant advantage of a rotating drum is that at low rota- 
tion rates it allows for the observation of many avalanches 
without the need to change the sample. After each 
avalanche the medium remains at rest relative to the 
drum while its surface angle is slowly increased by rota- 
tion, until it reaches ^max again. Then another avalanche 
occurs, and the process starts over. Thus it is possi- 
ble to record hundreds of avalanches, which is essential 
for performing a statistical analysis or for averaging out 
noise-like fluctuations in dynamical data. In addition to 
the avalanches observed at low rotations rates, the rotat- 
ing drum also allows the investigation of continuous flow 
down the slope at higher rotation rates. 

The sections of this paper cover different aspects of our 
rotating drum experiments. In Section we describe the 
experimental setup. Section addresses some issues ex- 
plored in previous studies of cohesive granular media in 
a rotating drum ^ ^ ^ , focussing on the surface 
angles of the medium before and after avalanches, and 
describing a statistical analysis of avalanche size based 
on these angles. In Section we investigate the transition 
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FIG. 1: Sketch of the experimental setup. 



between avalanches and continuous flow as a function of 
the liquid content, and map out the phase diagram of 
the system. Then a short section () is devoted to the 
morphology of the surface. In the rest of the paper we 
focus on characterizing the dynamics of cohesive flow. In 
Section we quantitatively investigate the flow dynamics 
during avalanches at different liquid contents by analyz- 
ing the time evolution of the averaged surface profile ob- 
tained from hundreds of avalanche events. In Section we 
analyze the dynamical properties of the continuous flow 
phase occurring at faster rotation. We pay special atten- 
tion to the nature of the viscoplastic flow, observed at 
the highest liquid contents which displays unique charac- 
teristics associated with coherent motion over the entire 
granular surface. In the final section wc summarize and 
discuss our results. 



DESCRIPTION OF THE APPARATUS 




FIG. 2: This is a reduced resolution figure Steps of im- 
age processing, (a) Snapshot taken without background il- 
lumination. The wet grains sticking to the glass walls make 
the surface profile hardly visible, (b) The same picture but 
with a circular lamp turned on in the background. In dif- 
fuse background illumination the monolayer of beads on the 
wall is much brighter than the 3.2 cm thick sample, (c) A 
brightness threshold is applied to the image to make the beads 
on the walls disappear, (d) Test of the algorithm: the sur- 
face data (black line) superimposed over a greyed version of 
the original photograph, (e) The surface profile with a fit- 
ted straight line determining the overall surface angle (9). 
(f) Pseudo-Sdimensional representation of the surface profile, 
shading indicates local angle (a(x)). This shading may be 
used to represent a given surface profile in a single line (see 
the horizontal stripe at the bottom). 



We studied glass spheres thoroughly mixed with small 
quantities of hydrocarbon oil. The viscosity of the oil 
used was 0.27 poise, and its surface tension was 0.02 N/m. 
The liquid content varied between r = 0.001 % and 5 % of 
the void volume. In this regime the flow of oil due to grav- 
ity can be neglected. Measurements were performed on 
three sizes of beads, with diameters d — 0.9 mm ± 11 %, 
d = 0.5 mm ±20% and d = 0.35 mm ±15%. Note the 
the beads differed not only in their size, but also pre- 
sumably in their microscopic surface structures, which 
are difficult to characterize or control. Nonetheless, most 
of the qualitative behavior is reproduced in the different 
grain samples, those differences which are signiflcant are 
mentioned below. 

The rotating drum (Fig. ^ was made of thick plexiglas, 
but the vertical walls were lined with glass plates in order 
to minimize electrostatic effects. To prevent slips along 
the circumference of our drum, we inserted a hollow, thin 
aluminum cylinder with a rough inner surface into the 
drum. The inner diameter of the drum was 16.8 cm, its 
width was 3.2 cm, and the granular filling was 30%. By 
performing measurements in a thinner (2 cm) drum we 
verified that, while wall effects are not negligible, they do 
not appear to modify the qualitative behavior. 

The drum was rotated by a computer-controlled pre- 
cision stepping motor which provided an extremely sta- 



ble rotation rate that could be varied in the range of 
~ 0.003 — 30 rotation per minute (rpm). The step size 
of rotation was 0.025°. The vibrations originating from 
the stepper motor were damped by multiple lead bricks 
attached to the apparatus. 

The experiment was recorded by a charge coupled de- 
vice video camera interfaced to a computer which could 
analyze the spatiotemporal evolution of the surface pro- 
file (height variations in the axial direction were negligi- 
ble). The steps of image processing are demonstrated in 
Fig. ||(a-c). When the grains are wet, the beads sticking 
to the glass walls make the surface hardly visible (Fig. 
H(a)). We overcame this problem by using background 
illumination (Fig. |^(b)) and applying a brightness thresh- 
old (Fig. 11(c)) on the image. We found the surface profile 
h{x) by an algorithm based on the detection of points of 
maximum contrast ]2^ . The resulting surface data are 
presented in Fig. ||(d-f ) in different ways. Fig. ^(d) shows 
a test of the algorithm: the detected surface is superim- 
posed on the photograph of the system as a thin black 
line. The limitations of our camera and computer gave us 
a resolution of 6t — 0.03 seconds and Sx = 0.5 mm. Fig. 
||(e) presents the h{x) curve together with a straight fine 
that was fitted to the curve in a rotation-invariant way: 
the squared sum of the geometrical distances of the sur- 
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face point from the line was minimized. The inchnation 
angle of this line gives the overall angle of the surface. 
This is the generalized version of the surface angle that 
has typically been used to describe the flat surface of 
dry samples. Finally, Fig. ||(f) shows the surface profile 
with a brightness-coded representation of the local angle 
{a{x) — arctan ) in the third dimension. The 

stripe at the bottom is a concise way of displaying this 
surface profile using the same shading. In sections and 
we use this representation of the surface topography to 
visualize the time evolution of the surface by placing the 
stripes corresponding to successive surface profiles under 
each other. 



DYNAMIC PROPERTIES FROM SURFACE 
ANGLE MEASUREMENTS 



The phenomenology of wet granular materials is much 
richer than that of dry samples. The attractive forces 
between the particles lead to correlated motion of the 
grains and qualitatively new avalanche dynamics. These 
new types of avalanches form rough and highly variable 
surfaces. In order to capture the most fundamental fea- 
tures of our system, however, we will initially neglect 
these surface features and examine only the overall sur- 
face angle (Fig. §(e)) to describe the observed surfaces 
with a single number. This approach allows us to com- 
pare our results directly to the behavior of dry materials, 
where the surface is virtually flat and is completely de- 
scribed by its inclination angle. 

Fig. H shows the time evolution of the surface angle 
for various rotation rates. Similarly to the dry materials 
we observe discrete avalanches at low rotation rates. We 
find that 6{t) in this regime is a typical sawtooth-shaped 
function looking very similar to the force-displacement 
curve of stick-slip processes, which is not very surprising 
since stick-slip-likc behavior is quite frequent in granular 
materials 
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. The local maximum and min- 
imum points of 9{t) correspond to the maximum angle 
of stability (^max) and the repose angle (6*^) respectively. 
Note, that if we plot the surface angle as a function of the 
rotation angle of the drum, then the increasing segments 
of the curve are straight lines with a slope of unity. If we 
gradually increase the angular velocity of the drum, then 
at a given rotation rate {fl = 4 rpm in Fig. ||) the flow 
becomes continuous, and after some transients the sur- 
face angle becomes constant. In this section we consider 
the avalanching regime (at = 0.12 rpm, except where 
noted otherwise) and investigate how the characteristic 
angles and the avalanche size distribution depends on the 
liquid content and the rotation rate. The transition to 
continuous flows is described in the subsequent section. 
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FIG. 3: This is a reduced resolution figure The temporal 
variation of the surface angle at various rotation rates (a) as 
a function of time, (b) as a function of the rotation angle of 
the drum. These data were taken using the large d = 0.9 mm 
beads, the liquid content was r = 0.009 %. Note that the 
curves are shifted by increments of 20° for clarity. 



The maximum and the repose angle 

Fig. ^ shows the maximum angle just before an 
avalanche, (^max), and the angle of repose after an 
avalanche, {9r), as a function of liquid content for the 
three different bead sizes (averaged over several hundred 
avalanches). The measured angles have distributions of 
finite width, the error bars in the figure indicate their 
standard deviation. These characteristic angles reveal 
essential properties of the material: ^max reflects the sta- 
bility of the sample while 9r is also related to the dis- 
sipation during the avalanche events. The various fea- 
tures of these curves reflect changes in the dynamical 
behavior corresponding to the three regimes mentioned 
in the introduction, which are described in this section 
through qualitative observations, and are more quantita- 
tively characterized in Section . 

The curves corresponding to the smaller bead sizes 
d = 0.5mm (Fig. |(b)) and d = 0.35mm (Fig. |(c)) are 
quite similar. At the lowest liquid contents, very small 
avalanches are observed. This is the granular regime, 
where the behavior is qualitatively similar to that of 
dry granular materials. Due to the interparticle cohe- 
sive forces, ^max in this regime is larger than the value 
measured in dry samples. On the other hand, 9r hardly 
changes, presumably because the presence of interstitial 
liquid does not significantly alter the dissipation in the 
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FIG. 4: The dependence of the maximum angle ^max and the 
repose angle 6,. on the liquid content for 3 different bead sizes, 
(a) d = 0.9 mm, SI = 0.12 rpm, (b) d = 0.5 mm, fl = 0.12 rpm, 
(c) d — 0.35 mm, Q — 0.03 rpm (a lower rotation rate was 
chosen, because at f2 = 0.12 rpm the smallest beads exhibit 
continuous flow for the lowest liquid contents). 



fluidized layer of rolling beads. 

At higher liquid contents, in the correlated regime, the 
increased cohesive force results in a larger increase in 
Smax- In this regime the beads on the top surface are 
more effectively stabilized by the cohesive forces, thus 
failure leading to an avalanche occurs in the bulk of the 
sample, and the avalanche flow includes the motion of 
correlated clumps. We find that 6',. also displays a sudden 
increase, but at somewhat higher liquid contents than 
the increase in 6',nax- This offset is possibly due to the 
motion of the clump fluidizing the underlying layer of 
beads, thus individual rolling grains help determine the 
final slope as in the granular regime. At higher liquid 
contents, the cohesive and viscous forces would prevent 
such fluidization leading to the increase in Or- 

At the highest liquid contents, 6'max slightly decreases 
in Fig. |^(a) and (b). This decrease in the stability of 
the material, also observed in our previous studies of the 
repose angle using a different technique is proba- 
bly related to lubrication effects decreasing the static 
friction among the grains. In this regime the nature 
of the avalanches changes again: the material is trans- 
ported by coherent flow extending over the whole sur- 
face of the sample. We did not observe this regime 
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FIG. 5: The average avalanche size and the width of the 
avalanche size distribution as a function of liquid content, 
calculated with two different methods, (a) Decrease in the 
surface angle: {A6), 2a (Ad), (b) Decrease in the potential 
energy: (AE), 2a{AE). The data were taken using medium 
size beads, d — 0.5 mm, Q = 0.12 rpm. 



for the d = 0.35 mm beads, possibly because, for those 
grains, the cohesive forces are much stronger relative to 
their weight. Careful examination of Fig. ^ reveals a 
few other differences between the different grains stud- 
ied such as the relatively large avalanches exhibited by 
the d = 0.9 mm beads even at very low liquid contents. 
These differences are difficult to interpret, since we only 
examined three grain samples which varied not only in 
their size, but presumably also in their microscopic sur- 
face properties. 



Statistics of avalanche sizes 

Thus far we have restricted our discussion to the aver- 
age values (6'max} and {9r)- However, the characteristic 
angles are infiuenced by the microstructure of the sam- 
ple, i.e., the actual configuration of the grains, which 
leads to inherent fluctuations and a non-trivial distribu- 
tion of the measured quantities. The analysis of these 
distributions may provide insight into the internal dy- 
namics of the material. In particular the statistics of 
avalanche sizes has attracted much recent attention due 
to its possible connection to self organized criticality 

0, 0, in, H, 11,0. 

In this section, we analyze the statistical distribution 
of avalanche sizes. The simplest quantity that describes 
the size of an avalanche is A9 = O^aax ~ which has 
been widely used to characterize avalanches among dry 

n|. In Fig. |(a) we plot 



grains in a rotating drum |17 
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FIG. 6: (color) Distribution of the avalanche size (AO) for 
various liquid contents, (a) Large beads, d — 0.9 mm, Q = 
0.12 rpm, (b) medium beads, d = 0.5 mm, Q = 0.12 rpm, (c) 
small beads, d = 0.35 mm, Q — 0.03 rpm. The righthand col- 
umn shows the number of avalanches included for each sam- 
ple - note that the statistics are much better for the large and 
medium beads than for the smaller beads. 



the average avalanche size 
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(1) 



(where A9i, i = L.A^ is the angle drop in the i 
avalanche) and the width of the avalanche size distri 
bution 



th 



N{N ~ 1) 



(2) 



for the d — 0.5 mm beads. (Note that we do not as- 
sume that the distribution is Gaussian, we just use the 
standard deviation as a measure of the width of the dis- 
tribution.) Since the definition of the surface angle is 
somewhat arbitrary in the case of our surfaces which 
can be both rough and have global curvature (as dis- 
cussed below), we also define an alternative measure of 
the avalanche size, AE, which equals the difference in 
the potential energy of the sample just before and just 
after the avalanche The analogous quantities (AE) 
and 2a{AE) are plotted in Fig. ||(b). The graphs indi- 
cate that the two different measures of the avalanche size 
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FIG. 7: Statistical measures of the avalanche size distribu- 
tions presented in Fig. pi for the 3 different bead sizes as 
indicated in the figure, (a-c) The average angle drop: {AO), 
and the standard deviation 2a{A6) (see text), (d-f) The most 
probable avalanche size, i.e. the maximum of P{A6) ("peak 
position"), and the width of the peaks (full width at half 
maximum, "F.W.H.M."). 



yield almost equivalent results, thus we use the standard 
A9 as the measure of the avalanche size in the follow- 
ing discussion. In order to characterize the distribution 
in more detail, in Fig. || we plot the measured proba- 
bility densities P{A9) for several liquid contents for the 
3 different bead sizes, where P{A9)d9 is the probability 
for the avalanche size to fall between A9 and A9 + d9 
(we estimate P{A9) from histograms). In Fig. |^ we also 
plot several statistical measures characterizing these dis- 
tributions including the peak position (the maximum of 
P{A9)) and the FWHM (full width of the peak at half 
of the maximum value). 

The statistical properties of the avalanches shown in 
Fig. H and Fig. |7| exhibit several interesting features. 
The size of the typical avalanche quickly increases with 
increasing liquid content reflecting the increasing cohe- 
sion, then a decrease is observable when the viscoplastic 
regime is reached. At the lowest liquid contents we ob- 
serve single peaks, i.e. the system has a characteristic 
avalanche size with noise-like variations. At larger liq- 
uid contents, however, in several cases we observe double 
peaks or extremely broad avalanche size distributions. 
This implies the existence of two or more characteristic 
avalanche sizes, i.e. both small and large avalanches may 
occur. A similar feature has been observed in fine cohe- 
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FIG. 8: Alternating small and large avalanches, (a) d = 
0.5mm, r = 0.125%, n = Q.4rpm, (b) d = 0.35mm, 
T = 0.062%, n = 0.15 rpm. (c) and (d) The initial (dashed 
line) and final (solid line) surface profiles for two successive 
avalanches. The parameters are the same as in (a) and (b) 
respectively. Unfortunately, due to the fiuctuations in the 
system these patterns of behavior were difficult to create re- 
producibly. 



sive powders ||2l|. As discussed below, this phenomenon 
is related to the fact that at high liquid contents the 
surface is far from being flat and both the limit of stabil- 
ity and the size of the avalanche depend strongly on the 
shape of the actual surface. 

Within a rotating drum experiment, the initial con- 
dition for an avalanche is the surface profile formed by 
the previous avalanche. The shapes of the surfaces after 
successive avalanches can thus be regarded as steps of a 
special dynamical map. In the granular regime the freely 
rolling grains always smooth out the surface, so this map- 
ping is trivial. However, in the correlated regime, where 
the existence of clumps leads to rugged surfaces, this map 
is rather complex. Due to the inherent fluctuations in the 
system, in most cases this mapping is stochastic, smaller 
and larger avalanches follow each other in a largely ran- 
dom manner. In some cases, however, the system reaches 
quasi-stable fixed points (when all the avalanches are very 
similar and their initial and final surfaces are rotated ver- 
sions of each other), or limit cycles, when the initial pro- 
file is recovered in two or more avalanches. Two examples 
for this latter case are shown in Fig. |[ Unfortunately the 
fluctuations in our system usually destroy regular behav- 
ior, and thus we could not perform a detailed investiga- 
tion of this map. But we expect that a corresponding 
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FIG. 9: (a-b) The dependence of the maximum angle Smax on 
the waiting time between successive avalanches. The system 
exhibits aging: longer waiting times between avalanches leads 
to increased stability. The continuous lines are logarithmic 
fits to the data, (a) d = 0.9 mm, (b) d = 0.5 mm. (c) The 
intensity of aging (i.e. the slope of the fitted straight lines 
in (a) and (b)) as a function of liquid content for the two 
different bead sizes. 



model system without strong fluctuations would display 
an extremely rich chaotic behavior. 



Logarithmic aging 

As expected from previous studies, and shown in Sec- 
tion , 6'niax strongly depends on the liquid content due 
to the increasing capillary forces. Our measurements 
indicate that 0max also depends on the rotation rate, 
il. Although (for small fi) the rotation rate only in- 
fluences the waiting time between successive avalanches, 
^max clearly decreases with increasing rotation rate. This 
effect (which is also visible in the measurements of the 
surface as a function of drum rotation rate shown be- 
low) is presented in Fig. ||(a) and (b), which show ^max 
as a function of the waiting time, i.e. the time between 
avalanche events which decreases with increasing rota- 
tion rate. The curves indicate a logarithmic increase in 
the stability of the system. In Fig. ^ (c) we also present 
the intensity of the aging process (quantifled as the slope 
of the fitted straight lines) . We do not observe aging for 
the dry samples, and the effect is somewhat stronger for 
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FIG. 10: This is a reduced resolution figure g^l- The oil at 
the contact point of a wet bead (similar to the ones used in 
the experiment, but d — 1.2 mm) and the glass plate on which 
it rests. The images are taken from above at different time 
instants via phase-contrast microscopy. The delay time from 
the instant when the bead was placed on the plate was (a) 
i Ri 5 s, (b) t ~ 15 s, (c) t « 25 s. The flow of oil to the contact 
point thus occurs on a timescale of several seconds, suggesting 
a natural explanation for our observed aging effects. 



the larger liquid contents, though on the basis of our data 
we cannot determine the exact functional form. 

Similar aging effects have been observed in several dif- 
ferent granular experiments ^, 33 1, but the underlying 
mechanisms are complex, since several factors may play 
an important role |3^ . Since in our experiment the effect 
is more pronounced at higher liquid contents, and the va- 
por pressure of oil is low, we attribute it to the motion 
of oil flowing towards the contact points rather than con- 
densation effects This assumption is also supported 
by our investigations of a wet contact by phase contrast 
microscopy. We placed a wet bead on the object plate 
and observed a flow of oil to the region of the contact 
point (Fig. |l^. Our observations revealed that a liq- 
uid bridge was formed instantaneously when the bead 
was placed on the plate, but the bridge's shape contin- 
ued changing slowly even after several seconds of waiting 
time. This finding further supports the importance of 
viscosity to the stability of wet granular media which 
was demonstrated by Samandani et al. ||, |^]. 



THE ROLE OF THE ROTATION RATE: 
AVALANCHES AND CONTINUOUS FLOW 

Hysteresis in the transition between continuous and 
avalanching flow 

As mentioned above, if the rotation rate is increased 
above a threshold the system exhibits continuous flow 
instead of discrete avalanches. Earlier studies in dry 
media revealed that the transition between avalanching 
and continuous flow is hysteretic in rotation rate ||l^, ^ . 
Now we have investigated the avalanching - continuous 
flow transition for various liquid contents. We slowly 
increased and then decreased the rotation rate and mea- 
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FIG. 11: (a) (b) (c) 6max and 9r as a function of rotation 
rate for 3 different liquid contents, d — 0.5 mm. The error 
bars represent the standard deviation of the observed angle 
distributions. 



sured the average surface angle in the continuous regime 
as well as ^max and 9r in the avalanching regime. The 
rotation rate was changed in discrete steps by factors 
between 1.2 and 1.4, each step was performed with a 
constant angular acceleration in 3 seconds. Between the 
steps the drum was rotated for typically 0.5 — 1 full ro- 
tation to determine whether the flow is continuous or 
avalanching (the minimum rotation was 90° at some of 
the lowest rotation rates, the maximum was 40 full rota- 
tions a,t — 30 rpm). We consider the flow continuous 
when the medium never comes to rest with respect to the 
drum, which coincides with A9 = 6'max ~ Or < 1°. The 
results of typical runs are presented in Fig. |ll|(a), (b) and 
(c), where ^max and Or are plotted as a function of the 
rotation rate. 

The three basic regimes of behavior observed pre- 
viously are also reflected in the nature of the transi- 
tion. In the granular regime at low liquid contents, (e.g. 
T = 0.037%), where the behavior is qualitatively similar 
to the dry case, we observe a clear hysteretic transition 
between continuous flow and avalanching (Fig. |ll|(a)). 
At a somewhat higher liquid content (r = 0.125%, see 
Fig. |ll|(b)) the correlated behavior is marked by a lack 
of hysteresis, and intermittent avalanches are observed 
at a relatively high rotation rate compared to the other 
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FIG. 12: (a) Phase diagram of the dynamical behavior as a 
function of the hquid content (r) and the rotation rate (0,). 
(a) large beads, d — 0.9 mm, (b) medium beads, d = 0.5 mm, 
(c) small beads, d = 0.35 mm. 



regimes. The continuous flow in this regime consists of 
a stream of separated chimps rather than the constant 
flux seen in the other two regimes. At r = 0.375 % (Fig. 
|TT|(c)) hysteresis is again observed, coinciding with the 
onset of the viscoplastic continuous flow, which is smooth 
and coherent over the entire sample. 

The origin of the hysteresis in the transition between 
continuous and avalanching flow is a fascinating subject 
0, m, H Hi, but a detailed discussion goes beyond 
the scope of this paper. We hypothesize that the hys- 
teresis is observed because the transitions on increasing 
and decreasing rotation rate have different origins, thus 
they do not necessarily coincide. A rough explanation of 



the transition from avalanching to continuous flow is that 
the successive avalanches merge. More precisely, at high 
velocities the rotation during the avalanche is becom- 
ing important and the constant energy supply prevents 
the grains from stopping. For the reverse transition, the 
continuous flow cannot be maintained at arbitrarily low 
velocities, because then the flowing layer cannot stay flu- 
idized, and the grains get stuck in small depressions of 
the underlying layer. 

The nature of the transition is naturally different in 
the different wetting regimes. In the granular and the 
viscoplastic regime the constant material flux in the con- 
tinuous flow phase is qualitatively different from the first 
increasing and then decreasing flux observed in the case 
of an avalanche. Thus the intermittent and the con- 
tinuous flow may represent two locally stable dynamical 
phases with a hysteretic transition between them. On the 
other hand in the correlated regime there is no qualitative 
change in the flow dynamics when the motion becomes 
continuous: the successive clumps do not merge even if 
their motion overlaps in time. According to the argu- 
ments advanced by Rajchenbach |jl^, the characteristic 
fall time of the grains is thus the same in the avalanch- 
ing and the continuous flow. Therefore we observe no 
hysteresis. 



The phase diagram 

We summarize our measurements of the different flow 
regimes in a phase diagram (Fig. |l^) indicating the na- 
ture of the flow as a function of liquid content and rota- 
tion rate. Each symbol corresponds to a measurement, 
grey symbols denote avalanches, black symbols denote 
continuous flow, white symbols mark the hysteretic re- 
gions. The most complete graph corresponds to the large, 
d — 0.9mm, beads (Fig. |l^(a)). Here the three regimes 
are clearly visible. Near r « 0.02 % the disappearance of 
hysteresis and the displacement of the transition point to- 
wards higher rotation rates indicates the transition from 
the granular to the correlated regime. The sudden change 
in the transition point around r « 0.1 % marks the onset 
of the viscoplastic continuous flow. 

The other two diagrams (Fig. |2|(b) and (c), medium 
{d — 0.5 mm) and small {d = 0.35 mm) beads, respec- 
tively) contain fewer points, but the main trends are vis- 
ible. A universal feature of the phase diagrams is that, 
upon approaching the correlated regime from the gran- 
ular one the transition point moves towards higher ro- 
tation rates. This suggests that the clump formation 
leads to more effective transport and quicker avalanches 
as was noted in our previous draining crater studies [ p^ . 
A general trend is that for smaller beads the transition 
to correlated behavior occurs at higher liquid contents, 
most probably due to the larger internal surface of these 
samples. 
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The most striking feature is that, as the bead size is 
decreased, the viscoplastic behavior and in particular the 
region of the viscoplastic continuous flow diminishes. For 
the d ~ 0.5 mm beads it is restricted to the region be- 
tween 0.3% < r < 0.8% and lrpm< Q < lOrpm, and 
the d = 0.35 mm beads exhibit no viscoplastic flow at all. 
A possible explanation is that, in samples consisting of 
smaller beads, the number of contacts at the boundary 
of a given volume of material is much larger, therefore 
the cohesive forces dominate inertial effects. For smaller 
beads the hysteresis of the continuous-to- avalanche tran- 
sition also seems to decrease, but this is difficult to inter- 
pret and will require further investigation to understand. 

We have performed some preliminary measurements 
on the role of the viscosity of the interstitial liquid. The 
results (not shown here) indicated that the fluid viscosity 
is indeed an important parameter, e.g. experiments with 
a highly viscous silicon oil yielded viscoplastic flow in a 
wider range of liquid contents and rotation rates. These 
results are also consistent with the findings of Samandani 
and KudroUi §. 



SURFACE MORPHOLOGY 

Up to this point we have described the surface by a 
single quantity, the overall surface angle, but the actual 
surfaces are rough and have significant curvature. This 
morphology contains additional information about the 
system, and should be useful in producing realistic mod- 
els of wet granular flow. Fig. |l^ shows sample surface 
profiles for various materials. The graph confirms that 
the wet samples form exotic surfaces, especially for the 
smaller beads. Note that, besides the overall shape of the 
surface, some smaller scale roughness is also observable 
for the larger liquid contents. These small scale features 
typically are random and fluctuate between avalanches. 
In the viscoplastic regime, however the coherence of the 
flow strongly reduces fluctuations, and a robust pattern 
is formed that is reproduced at the end of each avalanche 
(as discussed in section below). 

Since the surfaces of the largest {d — 0.9 mm) bead 
samples are relatively simple, we investigate them in 
more detail. Fig. |lj shows sample surfaces for the 
d = 0.9 mm beads at various liquid contents and ro- 
tation rates. In the granular regime the surface is al- 
most flat, correlated avalanches typically form concave 
surfaces, and the viscoplastic flow typically results in 
a convex surface. Furthermore, at the highest rotation 
rates, the surface is S-shaped due to inertial effects across 
the entire wetness range . The convexity of the sur- 
face can be characterized quantitatively several ways [ p7[ . 
A simple and robust measure is the integrated area (k) 
between the surface profile {h{x)) and the straight line 
connecting the two ends of the surface line. For a fiat 
surface where h(x) is a straight line, k = 0; k < im- 
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FIG. 13: This is a reduced resolution figure [^. Sample sur- 
faces for various materials, as indicated in the figure. The 
snapshots were taken right before an avalanche, the rotation 
rate was Q — 0.05 rpm for the large and medium beads and 
n — 0.043 rpm for the small beads. Note that the surface is 
close to flat in the granular regime, mostly concave in the cor- 
related regime and mostly convex in the viscoplastic regime. 
Note also the smaller scale structure at the smaller bead sizes. 
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FIG. 14; This is a reduced resolution figure Shape of 

the surface at various points of the parameter space for the 
large, d = 0.9 mm beads. A sketch of the phase diagram is 
shown in greyscale in the background: the dark grey region 
denotes continuous flow, the light grey region corresponds to 
avalanches; the white area is the hysteretic regime. The shape 
changes both as a function of Q and r. At the highest rotation 
rates the S-shaped regime ]17[ is clearly observable. 
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FIG. 15: The value of convexity parameter function of 

the rotation rate and the hquid content for the large beads, 
d = 0.9mm. Black symbols: concave surface (k < 0), white 
symbols: convex surface («: > 0), grey symbols: {k ~ 0). 
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FIG. 16: This is a reduced resolution figure Shape of 

the surface at various points of the parameter space, for the 
medium, d = 0.5 mm beads (see the caption of Fig. ^ for 
further explanation). The observed surface shapes are more 
complicated for these beads compared to the d = 0.9 mm 
beads, but the overall convex vs. concave nature is similar. 



plies a predominantly concave surface, k > implies a 
predominantly convex surface. 

In Fig. [ij the value of k is indicated by the color of 
symbols as a function of r and H.. Black symbols cor- 
respond to concave surfaces (negative k), and approxi- 
mately indicate the correlated regime. White symbols are 
used for positive k values, corresponding to the viscoplas- 



tic regime. Grey symbols indicate that k « 0, either 
because the surface is flat (upper left region of Fig. |l^ 
- this is approximately the granular regime), or because 
convex and concave parts are balanced (at the boundary 
of correlated and viscoplastic regimes). For comparison 
we also show the typical observed surface shapes for the 
medium size (d = 0.5mm) beads (Fig. |l^). Also here 
the convexity of the surface seems to correlate with the 
different flow regimes, however, some additional smaller 
scale features are also observable. 

Clearly, the phase boundaries in Fig. ^ are not ver- 
tical lines. This means that the liquid content is not 
the only parameter influencing the behavior: changing 
the rotation rate can also switch between different flow 
mechanisms. This observation implies that the formation 
of correlated clumps is a dynamical process with charac- 
teristic times comparable to the other timescales of the 
experiment. We expect that the different convexity in the 
different regimes could be reproduced relatively easily in 
the framework of a continuum model. A detailed expla- 
nation of the surface shapes, however, would probably 
require a full analysis of the avalanche dynamics. 



DETAILED DYNAMICS OF AVALANCHES 

The rotating drum apparatus allows us to obtain in- 
formation not just about the medium before and after 
the avalanche events, but also about the details of the 
grain motion during avalanche events. In order to ana- 
lyze the dynamics of avalanches we have obtained two- 
dimensional space-time matrices, h{x,t), characterizing 
the sample surfaces throughout the avalanche process 
which can then be analyzed to produce a variety of in- 
formation about the individual avalanches. By taking 
derivatives of the h(x,t) data, we obtain the local an- 
gle, a{x,t) = arctan[dxh{x,t)], and the local vertical ve- 
locity, u{x,t) = dth{x,t), of the surface profile. Fur- 
thermore, by then integrating the vertical velocity (using 
the continuity equation and assuming constant density 
psf), we also obtain the local flux in the avalanche, i.e. 
(t){x, t) — pw jl]j/2 9th{x' , t)dx' (where p is the grain den- 
sity, w is the width of the drum, and D is the drum di- 
ameter), which represents the material flowing through 
a given vertical plane at position x. 



17(a) 



In Fig. 

present snapshots of the progression of single avalanches 
for several typical liquid contents. Fig. |l^(b) displays 
a{x,t) as a function of space (horizontally) and time 
(downwards) for the same individual avalanches. In Fig. 
p7| (c-e), we present the average behavior of 300 — 500 
avalanches at the same liquid contents, and show sim- 
ilar graphs of the time evolution of {aix^t)), {u{x,t)), 
and {4>{x, t)), where () denotes averaging over avalanches. 
By obtaining these quantitative measures of the averaged 
properties, we can separate the robust characteristics of 
the avalanche dynamics from the large fluctuations which 
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FIG. 17: This is a reduced resolution figure [£2|. (color) 
Dynamics of avalanches of different types with grain size 
d = 0.5 mm, at = 0.12rpm. (a) Snapshots (at 0.1 — 0.2 s 
time intervals) of five single avalanches corresponding to dif- 
ferent liquid contents. The third dimension is used for a 
brightness-coded representation of the local slope (a). The 
green lines indicate approximate slip planes in the correlated 
regime, and the red line shows the traveling quasi-periodic 
surface features in the viscoplastic regime, (b) The local slope 
with the same brightness coding, as a function of space and 
time. A horizontal line corresponds to a surface profile at a 
given instant, and time increases downwards (thus avalanches 
propagate down and to the left). The slanted bright and dark 
regions correspond to the avalanche front and the kink respec- 
tively (see text). The stripes at higher liquid contents indicate 
lasting surface features, (c) (d) (e) The characteristic features 
of the avalanches averaged over 300 — 500 avalanche events. 
The displayed quantities are (c) the local surface angle, (d) 
the rate of change of local height, and (e) the local grain flux 
(the white lines indicate the point of maximum current), re- 
spectively. Note that the surface patterns at the highest oil 
contents are robust against averaging. 



are inherent in avalanche processes. 

The averaging process raises a few technical questions. 
First, the order of averaging and numerical differentia- 
tion is optional, i.e. we can either calculate a{x, t) and 
the other quantities from the raw h{x, t) data and per- 
form averaging afterwards, or take the averaged {h{x,t)) 
and calculate its numerical derivatives. We verified that 
the two methods yield equivalent results, and chose the 
former option. 



A second question is how to find the corresponding 
frames in different recorded avalanches. In order to make 
this question clear let us consider the averaging in more 
detail. Let us suppose that we have recorded the varia- 
tion of the surface profile for N avalanches, and we have 
N matrices hu'v^ i — 1..N denoting the local height in 
the ith avalanche at the position x = —D/2 + u5x and 
time t = Iq^ -\- V 5t, where 5x and 5t are the spatial and 
temporal resolutions respectively, and is "the starting 
time" of the ith avalanche. Then, to obtain the average 
behaviour, we calculate (/iu,d) = X^t^i hu]v/N . Since 
the properties of the avalanches (e.g. duration) vary be- 
tween avalanche events, the appropriate choice for t}^ is 
not trivial. We tried several algorithms and then decided 
to use the most robust one which was based on the vari- 
ation of the overall surface angle 9{t). We fitted straight 
lines to the segments of the sawtooth signal of 9{t) (see 
Fig. H), and calculated the time t^^^^ corresponding to the 
intersection of the straight lines at the beginning of an 



avalanche. Then we chose i, 



0.5 sec. 



A third problem arises from the fact that the hu',v ma- 
trices have empty elements: the surfaces do not extend 
from —D/2 to D/2. Simply excluding the empty ele- 
ments from the averaging leads to the appearance of some 
artifacts at the edges. To resolve this issue, we filled in 
the empty elements with the value of the closest elements 
in that row (as if the surface continued horizontally over 
the edge of the drum), performed the averaging, and then 
cut down those parts of the averaged profile that were 
outside the drum. Note that this procedure only influ- 
enced the behavior very close to the edges of the sample. 

In Fig. |l^we present the same data as series of curves, 
where the columns again represent different liquid con- 
tents. In the first row of graphs we plot the averaged 
surface profiles {{h(x))) at different time instants dur- 
ing the avalanche. The second row of graphs shows the 
variation of the local angle {{a{x))), and the third row 
plots the local flux ((^(z))) with triangles showing the 
position of maximum flux. In the following we analyze 
the dynamics in the 3 regimes based on Figs, a nd [l^ . 
We have also magnifled parts of Fig. |l^ in Fig. |ig| to em- 
phasize particular features in the data. For comparison 
we also show the most important avalanche types for the 
d — 0.9 mm beads in Fig. Apart from the absence of 
robust large-scale patterns for the d — 0.5 mm beads (dis- 
cussed below), the qualitative features of the avalanches 
are quite similar for the two types of beads although the 
effects are less pronounced for the d = 0.9 mm beads. 



The avalanche types 

Granular avalanches. With dry grains [ p^ p^ , p^ , |2g| , 
|30| , |35[ |, the surface remains almost flat throughout the 
avalanche, and the avalanches have a much longer dura- 
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FIG. 18: This is a reduced resolution figure [g2|. (color) Averaged dynamical parameters: {h{x)), {a{x)), and {cj>{x)) as a 
function of space at various time instants during the avalanche for 5 different liquid contents, d = 0.5 mm, = 0.12 rpm. The 
time interval between successive curves is 0.066 sec, some curves are marked with color as specified in the legend. In the lowest 
row, triangles mark the point of maximum current. The most important features for different oil contents are the following. 
Dry samples. The avalanches are very small, and our resolution is insufficient to resolve the detailed dynamics, t = 0.04 % 
Granular regime. The surface is close to flat, but the local angle is somewhat larger near the top of the slope. The dynamics is 
dominated by fronts of rolling grains. The region of steeper slope corresponding to the downward propagating front as well as 
the negative slope of the kink is visible in the graph of the local angle. The point of maximum current also moves downhill then 
uphill due to the fronts, r = 0.12 % Correlated regime: failure occurs along a slip plane and a block slides down. The block is 
best seen in the red curve, the edges are smoothed by the averaging over several hundred avalanches. Some rolling grains are 
still present, as marked by the presence of the kink, r — 0.37 % Correlated regime, usually with multiple avalanches, but these 
are averaged out. The fronts have disappeared, and the angle near the bottom of the slope decreases monotonically. There are 
lasting contacts during the avalanche, thus the grains cannot roll, and the material stops coherently as a block when it reaches 
the bottom, r = 5.00 % Viscoplastic regime. The surface moves coherently. The current extends over the whole surface during 
the whole duration of the avalanche. Curves of the local angle reveal the robust surface patterns as a series of local maxima. 
During the avalanche the local maxima move downwards (and spread out slightly), and a new maximum is formed near the 
top. 



tion and much smaller flux than in the wet media - as is 
expected due to the lack of cohesion. Our resolution is 
not sufficient to distinguish any propagating front in this 
case, but our system is very similar to that investigated 
by ||l2|, all the surface grains are close to the limit of 
their stability. Thus we assume that the failure mecha- 
nism is similar: the propagation of the front destabilizing 
the grains is quicker than the material flow. This type 
of front is marked only by a very slight dilation of the 
material, which we cannot detect. 

With the addition of very small amounts of liquid (e.g. 



T = 0.04%, see Fig. Fig. |l^(a)) the avalanches become 
much larger due to the onset of intergrain cohesion — 
this allows us to observe the dynamics in more detail. 
In this granular regime the avalanche is always initiated 
at the top of the surface, and the upper part of the sur- 
face is quickly destabilized. As the particles start moving 
downwards, a front of rolling grains travels downhill. The 
grains in the lower regions remain at rest until the rolling 
grains reach them. This type of front corresponds to the 
"start down" front of Douady et al. The difference 
in behavior compared to the perfectly dry material can 
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FIG. 19: This is a reduced resolution figure [g2[. (color) 
The various fronts observed in our setup in (a) the granu- 
lar regime (r = 0.04 %) and in (b) the low liquid content 
region of the correlated regime (r = 0.12 %), for d = 0.5 mm, 
Q = 0.12 rpm. The colors indicate the rate of change of local 
height during the avalanche (see Fig. ^) . We observe 4 differ- 
ent fronts. 1: destabilizing front. Propagates faster than the 
material flow via interactions between neighboring grains — 
destabilizes particles. 2: Avalanche front. Consists of rolling 
grains, which proceed on a medium at rest and destabilize 
particles by collisions. 3: Kink. The rolling grains accumu- 
late at the lower wall and form a region of low slope there. 
This slows down the incoming grains, thus this front travels 
uphill. 4: Rear front. The moving grains leave the static ones 
behind. 



probably be explained by the slight concave curvature of 
the surface: the grains lower on the slope are in a rela- 
tively stable state when the avalanche is initiated. 

When the downward-propagating front reaches the 
wall of the drum at the bottom of the slope, a kind of 
Shockwave is formed, and the rolling of grains is stopped 
by a region of smaller local angle travelling uphill. The 
second front is classified as a "stop up" front by , and 
corresponds to the "kink" seen in other experiments on 
dry beads (generally associated with individually rolling 
grains reaching a solid barrier at the bottom of the slope 

@ El)- 

Correlated Avalanches. At intermediate liquid con- 
tents, in the correlated regime (0.1% < t < 2%), 
the principal failure mechanism is a fracture along a 
curved slip plane (approximated by the green lines in Fig. 
|l^(a)), analogous to the dynamics of a class of geological 
events known as "slides" This is the failure mech- 

anism described by the Mohr-Coulomb model ||l|, |^. At 
T = 0.12% there is a single slip plane, but at larger liq- 
uid contents (r = 0.37%) the avalanches occur through 
a succession of local slip events. The medium becomes 



more cohesive with increasing liquid content, preventing 
grains from moving individually, and thus the kink dis- 
appears for T > 0.3% since the material moving as a 
connected block stops coherently when it hits the bot- 
tom. 

Viscoplastic Avalanches. The onset of the viscoplastic 
regime (r w 2 %) is accompanied by dramatic changes 
in the behavior. The flow becomes correlated across the 
entire granular surface as demonstrated by the parallel 
lines during the avalanche in Fig. ^ (b) and (b). Since 
the whole surface moves coherently (rather than break- 
ing apart and evolving separately in different parts of 
the surface), fluctuations are strongly suppressed [ p^ . 
This behavior is qualitatively similar to a different class 
of geological event, called a "debris flow" or "mudflow" 
p3| , PI The coherent nature of the motion in this 
regime leads to novel aspects of the dynamics which are 
not observed in other granular systems. 

Pattern formation in the viscoplastic regime 

One novel property of the viscoplastic avalanches is 
the robust topology of the top surface which sponta- 
neously forms a nearly periodic pattern (seen in Fig. 
|l^(a)). This surface structure is maintained essentially 
intact during the avalanche (note that the lines are con- 
tinuous throughout the avalanche in Fig. |l^(b)), indicat- 
ing that there are lasting contacts in the flowing layer. 
Moreover, the pattern is not random, but rather has fea- 
tures which are reproduced at the end of each avalanche. 
This is demonstrated most clearly in Fig. ^ and in Fig. 
|l^(b) , where the average of 347 avalanches of the r = 5 % 
sample has the same features as the typical individual 
avalanche shown in Fig. |l^(a) and (b). The robust na- 
ture of the surface structure of the wettest grains is in 
sharp contrast to the other regimes where averaging com- 
pletely smoothes out the smaller surface features. We 
can understand this behavior as resulting from coher- 
ence of the entire flow, which strongly reduces fluctu- 
ations in this regime. With minimal fluctuations, the 
final surface structure after each avalanche is essentially 
the same, thus setting the same initial condition for the 
next avalanche. With the same initial conditions for each 
avalanche, naturally the surface features are reproduced 
each time. 

Our experiments with the larger (d = 0.9 mm) beads 
also reveal some pattern formation (see Fig. ^(b)), but 
with a smaller characteristic size corresponding to 8 — 10 
grain diameters. The difference is probably due to the 
smaller ratio of the cohesive forces to the gravitational 
forces on the grains. Due to the small length scale of 
these patterns, they are largely averaged out, but we have 
also observed robust patterns that are repeated after each 
avalanche under several different sets of conditions (see 
e.g. Fig. |2^, T = 0.375%) (these other patterns are not 
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FIG. 20: This is a reduced resolution figure [£2|. (color) 
Avalanche dynamics in samples consisting of large beads 
d = 0.9 mm. (See also the caption of Fig. 



as periodic as at r = 5%, d = 0.5 mni). The dependence 
of the pattern on grain size and other factors will be the 
subject of a future investigation. 



DYNAMIC PROPERTIES OF CONTINUOUS 
FLOWS 

In the previous section we demonstrated that the 
avalanche dynamics are dramatically different in the 
three regimes of behavior. The different dynamics are 
also evident at higher rotation rates where the flow is 
continuous. In this section we compare the properties of 
the continuous flow in the three wetness regimes, with 
special attention to the viscoplastic flow. Here we only 
examine the d — 0.9 mm beads, since they exhibit vis- 
coplastic continuous flow in the widest range of parame- 
ter space. 

Fig. shows the evolution of the surface profile during 
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FIG. 21: This is a reduced resolution figure (color) Com- 
parison of the types of continuous flow at = 10 rpm (large 
beads: d = 0.9 mm), (a) snapshots at time intervals of 0.067s. 
The red marks indicate travelling surface features, (b) time 
evolution of the surface angle a{x,t). The travelling surface 
features appear as stripes. 



different types of continuous flow for 17 = 10 rpm. In the 
granular regime (e.g. r = and 0.002 %) the surface pro- 
file is practically a straight line with only noise- like varia- 
tions. Our resolution is insufficient to resolve the individ- 
ually rolling grains. By contrast, in the correlated regime 
the surface is clearly concave, and the ffow consists of a 
succession of clumps moving downhill. From the slope of 
the stripes in 21(b) we can estimate the velocity of the 
moving clumps: u « 50 — 60cm/s. In the viscoplastic 
regime the fiow is smooth again, and the travelling sur- 
face features demonstrate that there are lasting contacts 
in the flowing layer. The slopes of stripes indicate that 
the flow is much slower in this regime, u « 10 — 15 cm/s. 
We note that the rotation rate is the same for all the pre- 
sented liquid contents, meaning that the material flux is 
constant. The difference in surface velocity thus presum- 
ably indicates that the moving layer is much deeper in 
case of viscoplastic flow. 

In order to compare the viscoplastic flow to the or- 
dinary granular flow, we added some tracer particles to 
measure the surface velocity v more accurately as a func- 
tion of rotation rate, (Fig. ^). The viscoplastic flow 
is slower by a factor of two than the granular one, and 
more importantly, the two curves correspond to rather 
different functional forms. The curvature of the granu- 
lar curve suggest that at higher rotation rate the ffow 
depth is increased, as has been previously suggested p6| . 
Naturally in our finite setup the fiow depth cannot in- 
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FIG. 22: The velocity v at the surface as a function of ro- 
tation rate for the ordinary granular and the viscoplastic 
continuous flow (large beads, d = 0.9 mm, the liquid contents 
are r — 0.005 % and r = 0.25 % respectively). The continuous 
lines represent power-low fits to the data. The linear curve for 
the viscoplastic flow indicates that flow depth is independent 
of the rotation rate. 



crease infinitely which may explain why our measured 
v{Q) deviates from fi^/^ which would be expected from 
a linear velocity profile On the other hand, 

the recent work of Aranson et al. actually predicts a 2/3 
power law 
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The linear w(ri) for viscoplastic flow 
suggests a constant flow depth which is independent of 
the flow rate. This observation is consistent with the co- 
herent nature of this type of flow: the flow depth is not 
determined by local mechanisms |ll^ , but is fixed by the 
geometry of the whole system. In order to characterize 
the viscoplastic flow in greater detail and check whether 
the flow depth is really constant we have performed ex- 
plicit measurements to determine the extent and shape 
of the flowing layer. The results conflrmed that the flow 
depth is almost independent of the rotation rate p3|. 



CONCLUSIONS 

We have investigated many features of the dynamic 
behavior of wet granular samples in a rotating drum. We 
found that all these features are related and in particular 
all of them reflect the 3 fundamental regimes observed 
earlier |lO[ |. 

In the granular regime the cohesive forces are relatively 
small, and they have almost no effect in the fluidized 
flowing layer. Thus only a small increase is observable 
in ^max7 while 9r hardly changes. The grains are rolling 
freely, and avalanche dynamics is dominated by propa- 
gating fronts. These form nearly flat surfaces resulting 
in relatively small avalanche size fluctuations. At high ro- 
tation rates we observe continuous flow with a fluidized 
flowing layer, similarly to surface flows in dry samples. 

On the other hand, in the correlated regime the grains 
stick together and form clumps. The increased cohesion 



leads to a large increase in ^max and 9r ■ As the free rolling 
of grains becomes impossible, the avalanche fronts dis- 
appear and steep, rugged, and usually concave surfaces 
are formed. Variations of the surface shape lead to wide 
avalanche size distributions. In this regime the continu- 
ous flow consists of a succession of falling clumps, thus 
the hysteresis between avalanching and continuous flow 
disappears. 

Finally, in the viscoplastic regime the flow becomes co- 
herent over the entire sample, and 0max and A9 decrease 
slightly due to lubrication and viscous effects. The re- 
duced fluctuations lead to narrower avalanche size distri- 
butions and the formation of robust patterns that are re- 
produced at the end of each avalanche. The transition to 
continuous flow happens at relatively low rotation rates, 
and the hysteresis reappears. The coherence of the vis- 
coplastic continuous flow is manifested in traveling sur- 
face features and a velocity independent flow depth. 

While the flows we observe appear to have analogies 
with geological events, it is important to note that real 
geological materials usually consist of polydisperse irreg- 
ular particles often with very high (r « 100%) water 
contents and that the scaling of our system to geological 
lengthscales is non-trivial (indeed, based our on earlier 
work we expect that all of the results will be affected by 
the size of the container ||l^). Furthermore, avalanche 
studies in real soil have demonstrated additional phenom- 
ena associated with soil saturation |^ . Interestingly we 
still recover some of the basic dynamical processes in our 
model system, which should aid the description of qual- 
itatively different flow behaviors in the framework of a 
single model 0, ||. 

The changes in the dynamic behavior with wetting are 
associated with the increasingly coherent nature of the 
flow, i.e. the formation of coherently moving clusters - 
clumps - due to the increased cohesion and viscous ef- 
fects. Within a cluster, local velocity fluctuations should 
be suppressed, and thus the local granular temperature 
(T — (v^) — (v)^) should approach zero, but the clus- 
ters themselves both form and break apart during an 
avalanche process in a finite container. An important 
theoretical question raised by our data is how a length 
scale describing the size of the clumps may emerge from 
a granular flow model, and how such a length scale would 
vary with the type of media, the total size of the granular 
sample, the nature of intergranular adhesion, viscosity of 
the fluid, and the type of granular flow. While the present 
experiments have raised these questions, further experi- 
mental studies expanding the investigated phase space as 
well as computational modelling will almost certainly be 
necessary before a complete understanding is approach- 
able. 
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